Tracking arboviruses, their transmission vectors and potential hosts by nanopore sequencing of mosquitoes.

The risk to human health from mosquito-borne viruses such as dengue, chikungunya and yellow fever is increasing due to increased human expansion, deforestation and climate change. To anticipate and predict the spread and transmission of mosquito-borne viruses, a better understanding of the transmission cycle in mosquito populations is needed. We present a pathogen-agnostic combined sequencing protocol for identifying vectors, viral pathogens and their hosts or reservoirs using portable Oxford Nanopore sequencing. Using mosquitoes collected in São Paulo, Brazil, we extracted RNA for virus identification and DNA for blood meal and mosquito identification. Mosquitoes and blood meals were identified by comparing cytochrome c oxidase I (COI) sequences against a curated Barcode of Life Data System (BOLD). Viruses were identified using the SMART-9N protocol, which allows amplified DNA to be prepared with native barcoding for nanopore sequencing. Kraken 2 was employed to detect viral pathogens and Minimap2 and BOLD identified the contents of the blood meal. Due to the high similarity of some species, mosquito identification was conducted using blast after generation of consensus COI sequences using RACON polishing. This protocol can simultaneously uncover viral diversity, mosquito species and mosquito feeding habits. It also has the potential to increase understanding of mosquito genetic diversity and transmission dynamics of zoonotic mosquito-borne viruses.

Differences in the gut microbiota between Gurkhas and soldiers of British origin.

Previous work indicated that the incidence of travellers' diarrhoea (TD) is higher in soldiers of British origin, when compared to soldiers of Nepalese descent (Gurkhas). We hypothesise that the composition of the gut microbiota may be a contributing factor in the risk of developing TD in soldiers of British origin. This study aimed to characterise the gut microbial composition of Gurkha and non-Gurkha soldiers of the British Army. Recruitment of 38 soldiers (n = 22 Gurkhas, n = 16 non-Gurkhas) and subsequent stool collection, enabled shotgun metagenomic sequencing-based analysis of the gut microbiota. The microbiota of Gurkhas had significantly (P < 0.05) lower diversity, for both Shannon and Simpson diversity indices, using species level markers than the gut microbiota of non-Gurkha soldiers. Non-metric Multidimensional Scaling (NMDS) of the Bray-Curtis distance matrix revealed a significant difference in the composition of the gut microbiota between Gurkhas and non-Gurkha soldiers, at both the species level (P = 0.0178) and the genus level (P = 0.0483). We found three genera and eight species that were significantly enriched in the non-Gurkha group and one genus (Haemophilus) and one species (Haemophilus parainfluenzae) which were enriched in the Gurkha group. The difference in the microbiota composition between Gurkha soldiers and soldiers of British origin may contribute to higher colonization resistance against diarrhoeal pathogens in the former group. Our findings may enable further studies into interventions that modulate the gut microbiota of soldiers to prevent TD during deployment.

Genomic assessment of invasion dynamics of SARS-CoV-2 Omicron BA.1.

Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) variants of concern (VOCs) now arise in the context of heterogeneous human connectivity and population immunity. Through a large-scale phylodynamic analysis of 115,622 Omicron BA.1 genomes, we identified >6,000 introductions of the antigenically distinct VOC into England and analyzed their local transmission and dispersal history. We find that six of the eight largest English Omicron lineages were already transmitting when Omicron was first reported in southern Africa (22 November 2021). Multiple datasets show that importation of Omicron continued despite subsequent restrictions on travel from southern Africa as a result of export from well-connected secondary locations. Initiation and dispersal of Omicron transmission lineages in England was a two-stage process that can be explained by models of the country's human geography and hierarchical travel network. Our results enable a comparison of the processes that drive the invasion of Omicron and other VOCs across multiple spatial scales.

Multiplex MinION sequencing suggests enteric adenovirus F41 genetic diversity comparable to pre-COVID-19 era.

Human adenovirus F41 causes acute gastroenteritis in children, and has recently been associated with an apparent increase in paediatric hepatitis of unknown aetiology in the UK, with further cases reported in multiple countries. Relatively little is known about the genetic diversity of adenovirus F41 in UK children; and it is unclear what, if any, impact the COVID-19 pandemic has had on viral diversity in the UK. Methods that allow F41 to be sequenced from clinical samples without the need for viral culture are required to provide the genomic data to address these questions. Therefore, we evaluated an overlapping-amplicon method of sequencing adenovirus genomes from clinical samples using Oxford Nanopore technology. We applied this method to a small sample of adenovirus-species-F-positive extracts collected as part of standard care in the East of England region in January-May 2022. This method produced genomes with >75 % coverage in 13/22 samples and >50 % coverage in 19/22 samples. We identified two F41 lineages present in paediatric patients in the East of England in 2022. Where F41 genomes from paediatric hepatitis cases were available (n=2), these genomes fell within the diversity of F41 from the UK and continental Europe sequenced before and after the 2020-2021 phase of the COVID-19 pandemic. Our analyses suggest that overlapping amplicon sequencing is an appropriate method for generating F41 genomic data from high-virus-load clinical samples, and currently circulating F41 viral lineages were present in the UK and Europe before the COVID-19 pandemic.

Recurrent SARS-CoV-2 mutations in immunodeficient patients.

Long-term severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) infections in immunodeficient patients are an important source of variation for the virus but are understudied. Many case studies have been published which describe one or a small number of long-term infected individuals but no study has combined these sequences into a cohesive dataset. This work aims to rectify this and study the genomics of this patient group through a combination of literature searches as well as identifying new case series directly from the COVID-19 Genomics UK (COG-UK) dataset. The spike gene receptor-binding domain and N-terminal domain (NTD) were identified as mutation hotspots. Numerous mutations associated with variants of concern were observed to emerge recurrently. Additionally a mutation in the envelope gene, T30I was determined to be the second most frequent recurrently occurring mutation arising in persistent infections. A high proportion of recurrent mutations in immunodeficient individuals are associated with ACE2 affinity, immune escape, or viral packaging optimisation. There is an apparent selective pressure for mutations that aid cell-cell transmission within the host or persistence which are often different from mutations that aid inter-host transmission, although the fact that multiple recurrent de novo mutations are considered defining for variants of concern strongly indicates that this potential source of novel variants should not be discounted.

Context-specific emergence and growth of the SARS-CoV-2 Delta variant.

The SARS-CoV-2 Delta (Pango lineage B.1.617.2) variant of concern spread globally, causing resurgences of COVID-19 worldwide1,2. The emergence of the Delta variant in the UK occurred on the background of a heterogeneous landscape of immunity and relaxation of non-pharmaceutical interventions. Here we analyse 52,992 SARS-CoV-2 genomes from England together with 93,649 genomes from the rest of the world to reconstruct the emergence of Delta and quantify its introduction to and regional dissemination across England in the context of changing travel and social restrictions. Using analysis of human movement, contact tracing and virus genomic data, we find that the geographic focus of the expansion of Delta shifted from India to a more global pattern in early May 2021. In England, Delta lineages were introduced more than 1,000 times and spread nationally as non-pharmaceutical interventions were relaxed. We find that hotel quarantine for travellers reduced onward transmission from importations; however, the transmission chains that later dominated the Delta wave in England were seeded before travel restrictions were introduced. Increasing inter-regional travel within England drove the nationwide dissemination of Delta, with some cities receiving more than 2,000 observable lineage introductions from elsewhere. Subsequently, increased levels of local population mixing-and not the number of importations-were associated with the faster relative spread of Delta. The invasion dynamics of Delta depended on spatial heterogeneity in contact patterns, and our findings will inform optimal spatial interventions to reduce the transmission of current and future variants of concern, such as Omicron (Pango lineage B.1.1.529).

Molecular and genomic investigation of an urban outbreak of dengue virus serotype 2 in Angola, 2017-2019.

BACKGROUND: The transmission patterns and genetic diversity of dengue virus (DENV) circulating in Africa remain poorly understood. Circulation of the DENV serotype 1 (DENV1) in Angola was detected in 2013, while DENV serotype 2 (DENV2) was detected in 2018. Here, we report results from molecular and genomic investigations conducted at the Ministry of Health national reference laboratory (INIS) in Angola on suspected dengue cases detected between January 2017 and February 2019. METHODS: A total of 401 serum samples from dengue suspected cases were collected in 13 of the 18 provinces in Angola. Of those, 351 samples had complete data for demographic and epidemiological analysis, including age, gender, province, type of residence, clinical symptoms, as well as dates of onset of symptoms and sample collection. RNA was extracted from residual samples and tested for DENV-RNA using two distinct real time RT-PCR protocols. On-site whole genome nanopore sequencing was performed on RT-PCR+ samples. Bayesian coalescent models were used to estimate date and origin of outbreak emergence, as well as population growth rates. RESULTS: Molecular screening showed that 66 out of 351 (19%) suspected cases were DENV-RNA positive across 5 provinces in Angola. DENV RT-PCR+ cases were detected more frequently in urban sites compared to rural sites. Of the DENV RT-PCR+ cases most were collected within 6 days of symptom onset. 93% of infections were confirmed by serotype-specific RT-PCR as DENV2 and 1 case (1.4%) was confirmed as DENV1. Six CHIKV RT-PCR+ cases were also detected during the study period, including 1 co-infection of CHIKV with DENV1. Most cases (87%) were detected in Luanda during the rainy season between April and October. Symptoms associated with severe dengue were observed in 11 patients, including 2 with a fatal outcome. On-site nanopore genome sequencing followed by genetic analysis revealed an introduction of DENV2 Cosmopolitan genotype (also known as DENV2-II genotype) possibly from India in or around October 2015, at least 1 year before its detection in the country. Coalescent models suggest relatively moderately rapid epidemic growth rates and doubling times, and a moderate expansion of DENV2 in Angola during the studied period. CONCLUSION: This study describes genomic, epidemiological and demographic characteristic of predominately urban transmission of DENV2 in Angola. We also find co-circulation of DENV2 with DENV1 and CHIKV and report several RT-PCR confirmed severe dengue cases in the country. Increasing dengue awareness in healthcare professional, expanding the monitorization of arboviral epidemics across the country, identifying most common mosquito breeding sites in urban settings, implementing innovative vector control interventions and dengue vaccination campaigns could help to reduce vector presence and DENV transmission in Angola.

Generation and transmission of interlineage recombinants in the SARS-CoV-2 pandemic.

We present evidence for multiple independent origins of recombinant SARS-CoV-2 viruses sampled from late 2020 and early 2021 in the United Kingdom. Their genomes carry single-nucleotide polymorphisms and deletions that are characteristic of the B.1.1.7 variant of concern but lack the full complement of lineage-defining mutations. Instead, the remainder of their genomes share contiguous genetic variation with non-B.1.1.7 viruses circulating in the same geographic area at the same time as the recombinants. In four instances, there was evidence for onward transmission of a recombinant-origin virus, including one transmission cluster of 45 sequenced cases over the course of 2 months. The inferred genomic locations of recombination breakpoints suggest that every community-transmitted recombinant virus inherited its spike region from a B.1.1.7 parental virus, consistent with a transmission advantage for B.1.1.7's set of mutations.

A multi-targeting pre-clinical candidate against drug-resistant tuberculosis.

FNDR-20081 [4-{4-[5-(4-Isopropyl-phenyl)- [1,2,4]oxadiazol-3-ylmethyl]-piperazin-1-yl}-7-pyridin-3-yl-quinoline] is a novel, first in class anti-tubercular pre-clinical candidate against sensitive and drug-resistant Mycobacterium tuberculosis (Mtb). In-vitro combination studies of FNDR-20081 with first- and second-line drugs exhibited no antagonism, suggesting its compatibility for developing new combination-regimens. FNDR-20081, which is non-toxic with no CYP3A4 liability, demonstrated exposure-dependent killing of replicating-Mtb, as well as the non-replicating-Mtb, and efficacy in a mouse model of infection. Whole genome sequencing (WGS) of FNDR-20081 resistant mutants revealed the identification of pleotropic targets: marR (Rv0678), a regulator of MmpL5, a transporter/efflux pump mechanism for drug resistance; and Rv3683, a putative metalloprotease potentially involved in peptidoglycan biosynthesis. In summary, FNDR-20081 is a promising first in class compound with the potential to form a new combination regimen for MDR-TB treatment.

CLIMB-COVID: continuous integration supporting decentralised sequencing for SARS-CoV-2 genomic surveillance.

In response to the ongoing SARS-CoV-2 pandemic in the UK, the COVID-19 Genomics UK (COG-UK) consortium was formed to rapidly sequence SARS-CoV-2 genomes as part of a national-scale genomic surveillance strategy. The network consists of universities, academic institutes, regional sequencing centres and the four UK Public Health Agencies. We describe the development and deployment of CLIMB-COVID, an encompassing digital infrastructure to address the challenge of collecting and integrating both genomic sequencing data and sample-associated metadata produced across the COG-UK network.

Spatiotemporal invasion dynamics of SARS-CoV-2 lineage B.1.1.7 emergence.

Understanding the causes and consequences of the emergence of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) variants of concern is crucial to pandemic control yet difficult to achieve because they arise in the context of variable human behavior and immunity. We investigated the spatial invasion dynamics of lineage B.1.1.7 by jointly analyzing UK human mobility, virus genomes, and community-based polymerase chain reaction data. We identified a multistage spatial invasion process in which early B.1.1.7 growth rates were associated with mobility and asymmetric lineage export from a dominant source location, enhancing the effects of B.1.1.7's increased intrinsic transmissibility. We further explored how B.1.1.7 spread was shaped by nonpharmaceutical interventions and spatial variation in previous attack rates. Our findings show that careful accounting of the behavioral and epidemiological context within which variants of concern emerge is necessary to interpret correctly their observed relative growth rates.

Assessing transmissibility of SARS-CoV-2 lineage B.1.1.7 in England.

The SARS-CoV-2 lineage B.1.1.7, designated variant of concern (VOC) 202012/01 by Public Health England1, was first identified in the UK in late summer to early autumn 20202. Whole-genome SARS-CoV-2 sequence data collected from community-based diagnostic testing for COVID-19 show an extremely rapid expansion of the B.1.1.7 lineage during autumn 2020, suggesting that it has a selective advantage. Here we show that changes in VOC frequency inferred from genetic data correspond closely to changes inferred by S gene target failures (SGTF) in community-based diagnostic PCR testing. Analysis of trends in SGTF and non-SGTF case numbers in local areas across England shows that B.1.1.7 has higher transmissibility than non-VOC lineages, even if it has a different latent period or generation time. The SGTF data indicate a transient shift in the age composition of reported cases, with cases of B.1.1.7 including a larger share of under 20-year-olds than non-VOC cases. We estimated time-varying reproduction numbers for B.1.1.7 and co-circulating lineages using SGTF and genomic data. The best-supported models did not indicate a substantial difference in VOC transmissibility among different age groups, but all analyses agreed that B.1.1.7 has a substantial transmission advantage over other lineages, with a 50% to 100% higher reproduction number.

Transmission, adaptation and geographical spread of the Pseudomonas aeruginosa Liverpool epidemic strain.

The Liverpool epidemic strain (LES) is an important transmissible clonal lineage of Pseudomonas aeruginosa that chronically infects the lungs of people with cystic fibrosis (CF). Previous studies have focused on the genomics of the LES in a limited number of isolates, mostly from one CF centre in the UK, and from studies highlighting identification of the LES in Canada. Here we significantly extend the current LES genome database by genome sequencing 91 isolates from multiple CF centres across the UK, and we describe the comparative genomics of this large collection of LES isolates from the UK and Canada. Phylogenetic analysis revealed that the 145 LES genomes analysed formed a distinct clonal lineage when compared with the wider P. aeruginosa population. Notably, the isolates formed two clades: one associated with isolates from Canada, and the other associated with UK isolates. Further analysis of the UK LES isolates revealed clustering by clinic geography. Where isolates clustered closely together, the association was often supported by clinical data linking isolates or patients. When compared with the earliest known isolate, LESB58 (from 1988), many UK LES isolates shared common loss-of-function mutations, such as in genes gltR and fleR. Other loss-of-function mutations identified in previous studies as common adaptations during CF chronic lung infections were also identified in multiple LES isolates. Analysis of the LES accessory genome (including genomic islands and prophages) revealed variations in the carriage of large genomic regions, with some evidence for shared genomic island/prophage complement according to clinic location. Our study reveals divergence and adaptation during the spread of the LES, within the UK and between continents.

Evaluation of full-length nanopore 16S sequencing for detection of pathogens in microbial keratitis.

BACKGROUND: Microbial keratitis is a leading cause of preventable blindness worldwide. Conventional sampling and culture techniques are time-consuming, with over 40% of cases being culture-negative. Nanopore sequencing technology is portable and capable of generating long sequencing reads in real-time. The aim of this study is to evaluate the potential of nanopore sequencing directly from clinical samples for the diagnosis of bacterial microbial keratitis. METHODS: Using full-length 16S rRNA amplicon sequences from a defined mock microbial community, we evaluated and benchmarked our bioinformatics analysis pipeline for taxonomic assignment on three different 16S rRNA databases (NCBI 16S RefSeq, RDP and SILVA) with clustering at 97%, 99% and 100% similarities. Next, we optimised the sample collection using an ex vivo porcine model of microbial keratitis to compare DNA recovery rates of 12 different collection methods: 21-gauge needle, PTFE membrane (4 mm and 6 mm), Isohelix™ SK-2S, Sugi® Eyespear, Cotton, Rayon, Dryswab™, Hydraflock®, Albumin-coated, Purflock®, Purfoam and Polyester swabs. As a proof-of-concept study, we then used the sampling technique that provided the highest DNA recovery, along with the optimised bioinformatics pipeline, to prospectively collected samples from patients with suspected microbial keratitis. The resulting nanopore sequencing results were then compared to standard microbiology culture methods. RESULTS: We found that applying alignment filtering to nanopore sequencing reads and aligning to the NCBI 16S RefSeq database at 100% similarity provided the most accurate bacterial taxa assignment. DNA concentration recovery rates differed significantly between the collection methods (p < 0.001), with the Sugi® Eyespear swab providing the highest mean rank of DNA concentration. Then, applying the optimised collection method and bioinformatics pipeline directly to samples from two patients with suspected microbial keratitis, sequencing results from Patient A were in agreement with culture results, whilst Patient B, with negative culture results and previous antibiotic use, showed agreement between nanopore and Illumina Miseq sequencing results. CONCLUSION: We have optimised collection methods and demonstrated a novel workflow for identification of bacterial microbial keratitis using full-length 16S nanopore sequencing.

Establishment and lineage dynamics of the SARS-CoV-2 epidemic in the UK.

The United Kingdom's COVID-19 epidemic during early 2020 was one of world's largest and was unusually well represented by virus genomic sampling. We determined the fine-scale genetic lineage structure of this epidemic through analysis of 50,887 severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) genomes, including 26,181 from the UK sampled throughout the country's first wave of infection. Using large-scale phylogenetic analyses combined with epidemiological and travel data, we quantified the size, spatiotemporal origins, and persistence of genetically distinct UK transmission lineages. Rapid fluctuations in virus importation rates resulted in >1000 lineages; those introduced prior to national lockdown tended to be larger and more dispersed. Lineage importation and regional lineage diversity declined after lockdown, whereas lineage elimination was size-dependent. We discuss the implications of our genetic perspective on transmission dynamics for COVID-19 epidemiology and control.

Context-specific emergence and growth of the SARS-CoV-2 Delta variant.

The Delta variant of concern of SARS-CoV-2 has spread globally causing large outbreaks and resurgences of COVID-19 cases. The emergence of Delta in the UK occurred on the background of a heterogeneous landscape of immunity and relaxation of non-pharmaceutical interventions. Here we analyse 52,992 Delta genomes from England in combination with 93,649 global genomes to reconstruct the emergence of Delta, and quantify its introduction to and regional dissemination across England, in the context of changing travel and social restrictions. Through analysis of human movement, contact tracing, and virus genomic data, we find that the focus of geographic expansion of Delta shifted from India to a more global pattern in early May 2021. In England, Delta lineages were introduced >1,000 times and spread nationally as non-pharmaceutical interventions were relaxed. We find that hotel quarantine for travellers from India reduced onward transmission from importations; however the transmission chains that later dominated the Delta wave in England had been already seeded before restrictions were introduced. In England, increasing inter-regional travel drove Delta's nationwide dissemination, with some cities receiving >2,000 observable lineage introductions from other regions. Subsequently, increased levels of local population mixing, not the number of importations, was associated with faster relative growth of Delta. Among US states, we find that regions that previously experienced large waves also had faster Delta growth rates, and a model including interactions between immunity and human behaviour could accurately predict the rise of Delta there. Delta’s invasion dynamics depended on fine scale spatial heterogeneity in immunity and contact patterns and our findings will inform optimal spatial interventions to reduce transmission of current and future VOCs such as Omicron.

Context-specific emergence and growth of the SARS-CoV-2 Delta variant.

The Delta variant of concern of SARS-CoV-2 has spread globally causing large outbreaks and resurgences of COVID-19 cases 1-3 . The emergence of Delta in the UK occurred on the background of a heterogeneous landscape of immunity and relaxation of non-pharmaceutical interventions 4,5 . Here we analyse 52,992 Delta genomes from England in combination with 93,649 global genomes to reconstruct the emergence of Delta, and quantify its introduction to and regional dissemination across England, in the context of changing travel and social restrictions. Through analysis of human movement, contact tracing, and virus genomic data, we find that the focus of geographic expansion of Delta shifted from India to a more global pattern in early May 2021. In England, Delta lineages were introduced >1,000 times and spread nationally as non-pharmaceutical interventions were relaxed. We find that hotel quarantine for travellers from India reduced onward transmission from importations; however the transmission chains that later dominated the Delta wave in England had been already seeded before restrictions were introduced. In England, increasing inter-regional travel drove Delta's nationwide dissemination, with some cities receiving >2,000 observable lineage introductions from other regions. Subsequently, increased levels of local population mixing, not the number of importations, was associated with faster relative growth of Delta. Among US states, we find that regions that previously experienced large waves also had faster Delta growth rates, and a model including interactions between immunity and human behaviour could accurately predict the rise of Delta there. Delta's invasion dynamics depended on fine scale spatial heterogeneity in immunity and contact patterns and our findings will inform optimal spatial interventions to reduce transmission of current and future VOCs such as Omicron.

Rapid viral metagenomics using SMART-9N amplification and nanopore sequencing.

Emerging and re-emerging viruses are a global health concern. Genome sequencing as an approach for monitoring circulating viruses is currently hampered by complex and expensive methods. Untargeted, metagenomic nanopore sequencing can provide genomic information to identify pathogens, prepare for or even prevent outbreaks. SMART (Switching Mechanism at the 5' end of RNA Template) is a popular approach for RNA-Seq but most current methods rely on oligo-dT priming to target polyadenylated mRNA molecules. We have developed two random primed SMART-Seq approaches, a sequencing agnostic approach 'SMART-9N' and a version compatible rapid adapters  available from Oxford Nanopore Technologies 'Rapid SMART-9N'. The methods were developed using viral isolates, clinical samples, and compared to a gold-standard amplicon-based method. From a Zika virus isolate the SMART-9N approach recovered 10kb of the 10.8kb RNA genome in a single nanopore read. We also obtained full genome coverage at a high depth coverage using the Rapid SMART-9N, which takes only 10 minutes and costs up to 45% less than other methods. We found the limits of detection of these methods to be 6 focus forming units (FFU)/mL with 99.02% and 87.58% genome coverage for SMART-9N and Rapid SMART-9N respectively. Yellow fever virus plasma samples and SARS-CoV-2 nasopharyngeal samples previously confirmed by RT-qPCR with a broad range of Ct-values were selected for validation. Both methods produced greater genome coverage when compared to the multiplex PCR approach and we obtained the longest single read of this study (18.5 kb) with a SARS-CoV-2 clinical sample, 60% of the virus genome using the Rapid SMART-9N method. This work demonstrates that SMART-9N and Rapid SMART-9N are sensitive, low input, and long-read compatible alternatives for RNA virus detection and genome sequencing and Rapid SMART-9N improves the cost, time, and complexity of laboratory work.

Evaluating the Effects of SARS-CoV-2 Spike Mutation D614G on Transmissibility and Pathogenicity.

Global dispersal and increasing frequency of the SARS-CoV-2 spike protein variant D614G are suggestive of a selective advantage but may also be due to a random founder effect. We investigate the hypothesis for positive selection of spike D614G in the United Kingdom using more than 25,000 whole genome SARS-CoV-2 sequences. Despite the availability of a large dataset, well represented by both spike 614 variants, not all approaches showed a conclusive signal of positive selection. Population genetic analysis indicates that 614G increases in frequency relative to 614D in a manner consistent with a selective advantage. We do not find any indication that patients infected with the spike 614G variant have higher COVID-19 mortality or clinical severity, but 614G is associated with higher viral load and younger age of patients. Significant differences in growth and size of 614G phylogenetic clusters indicate a need for continued study of this variant.

Improvements to the ARTIC multiplex PCR method for SARS-CoV-2 genome sequencing using nanopore.

Genome sequencing has been widely deployed to study the evolution of SARS-CoV-2 with more than 90,000 genome sequences uploaded to the GISAID database. We published a method for SARS-CoV-2 genome sequencing (https://www.protocols.io/view/ncov-2019-sequencing-protocol-bbmuik6w) online on January 22, 2020. This approach has rapidly become the most popular method for sequencing SARS-CoV-2 due to its simplicity and cost-effectiveness. Here we present improvements to the original protocol: i) an updated primer scheme with 22 additional primers to improve genome coverage, ii) a streamlined library preparation workflow which improves demultiplexing rate for up to 96 samples and reduces hands-on time by several hours and iii) cost savings which bring the reagent cost down to £10 per sample making it practical for individual labs to sequence thousands of SARS-CoV-2 genomes to support national and international genomic epidemiology efforts.